Gravitational waves from binaries on unbound orbits 
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A generalized true anomaly-type parametrization, convenient to describe both bound and open 
orbits of a two-body system in general relativity is introduced. A complete description of the time 
evolution of both the radial and of the angular equations of a binary system taking into account the 
first order post-newtonian (1PN) is given. The gravitational radiation field emitted by the system 
is computed in the 1PN approximation including higher multipole moments beyond the standard 
quadrupole term. The gravitational waveforms in the time domain are explicitly given up to the 
1PN order for unbound orbits, but the results are also illustrated on binaries on elliptic orbits with 
special attention given to the effects of eccentricity. 

PACS numbers: 



I. INTRODUCTION 



o 

Oh 

<D 

Among the significant sources of gravitational waves (GWs) for ground-based interferometers [lHJ] are parabolic and 
hyperbolic encounters of compact objects. Similarly to eccentric binaries these sources have an analytic description 
which allows one to evaluate the emitted gravitational wave signal in the post- Newtonian (PN) approximation [f|. 
During the close encounter of black holes sufficient amount of energy can be radiated in GWs resulting in the formation 
of a binary system. These phenomena are analyzed in Q with the method of perturbation of orbital elements. Close 
encounter type sources provide more characteristic wave signals than typical elliptic binaries, since after the first 
quasiparabolic burstlike signal it emits further bursts quasiperiodically, until it reaches the validity region of the 
elliptic description. Recent estimates show that these binary interactions can contribute significantly to the sources 
of forthcoming advanced gravitational wave observatories |7j . 

A new generation of ground-based interferometric detectors and the Laser Interferometer Space Antenna (LISA) 
increase the possibility of direct detection of GWs. The LISA project aims to install a detector system into space 
and extend the GW observations to the low-frequency range which is inaccessible to ground-based interferometers due 
to the local gravitational noise. Depending on the system parameters and initial conditions parabolic and hyperbolic 
i encounters produce GWs in the sensitivity range of LISA, between 10 -4 and 10 _1 Hz, as it is indicated in our study 
of specific binary systems. There are two well-known sources, the Hulse- Taylor pulsar [l(| and the first double pulsar 
system J0737-3039 [ll|, [HJ representing unique laboratories for relativistic gravitational physics which provide GWs 
in the frequency range of LISA. The computations presented here can be applied to the advanced ground-based 
detectors, too. 

The classical motion and the encounter of compact objects are usually considered in the weak field, post-Newtonian 
regime, where the motion of the binary system is well described by perturbed Keplerian orbits. The equations 
of motion and the emitted radiation of these systems are analyzed in detail in the literature. First order (1PN) 
contributions to the dynamical evolution for general, eccentric orbits are discussed in [l3j . Moreover, the solution 
for the radial and angular motion in quasi-Newtonian parametric form for closed and open orbits was also presented. 
The contributions of 2PN effects to the periastron advance and the orbital period are given in [l4j], where spin-orbit 
contribution to the secular precession of the orbit is also analyzed. Spin and tail effects on the energy and angular 
momentum losses were determined in (l5| . The energy and angular momentum fluxes, moreover the evolution of the 
orbital elements and the emitted gravitational waveform of compact binaries to 2PN order are described by [l6[ for 
eccentric orbits. Their results are expressed by the generalized quasi-Keplerian description and the eccentric anomaly 
parametrization for general elliptic binary orbits. Later, the PN description of binaries was extended to higher orders 
and the energy and angular momentum fluxes are given for inspiralling compact binaries in quasielliptical orbits at 
3PN [13, Ell- Moreover, 3.5PN corrections to the acceleration for a spherically symmetric nonrotating self-gravitating 
star can be found, e.g. in (l9| . The 1PN accurate frequency domain waveforms are given recently in terms of tensor 
spherical harmonics for eccentric nonspinning compact binaries in [20J. 

The main expressions for the Keplerian motion of binaries on unbound orbits have been derived by many authors in 
a variety of different contexts. The energy and angular momentum losses for hyperbolic encounters and the capture 
of particles arriving from infinity were analyzed by Hansen [2l| , whose work is an extension of the classical results 
of Peters and Mathews [22j for elliptic orbits. In [23] the gravitational waveforms and the energy spectrum from 
two point masses in arbitrary unbound orbits were given in the quadrupole formalism. As a result, the waveform is 
found to be highly peaked near the periastron for large eccentricities. Gravitational bremsstrahlung radiation has 
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been studied in pij by considering stars interacting on unbound orbits for the case of large impact parameters in the 
postlinear formalism. This description is valid for arbitrary stellar velocities and mass ratio of the two stars. Recently, 
the strain amplitude and the radiation luminosity have been analyzed in [25j for binaries where the periastron distance 
is much larger than the Schwarzschild radius of the individual stars. 

There are detailed studies of the dynamics of black holes in a stellar cluster using numerical simulations that include 
the effects of gravitational radiation. Cross sections for mergers resulting from gravitational radiation during two and 
three-body encounters for a range of binary semimajor axes and mass ratios were presented in (26| with the use of 
numerical techniques. In [27] fluxes of energy and angular momentum and the gravitational waveforms produced by 
a point particle traveling around a Schwarzschild black hole on arbitrary bound and unbound orbits are given. 

In our work we consider binary systems of compact objects in which the components are assumed to have negligible 
rotation. Both the classical motion (i.e. ignoring radiation reaction) and the radiation field itself are described in the 
1PN approximation including higher multipole moments beyond the quadrupole term. We introduce a generalized 
true anomaly parametrization to describe the radial motion and to express all relevant dynamical quantities in terms 
of it. This new definition of the radial parametrization makes a unified description of both open and closed orbits 
possible, i.e. the parallel discussion of the elliptic, parabolic and hyperbolic cases. Moreover, we integrate the time 
dependence of the generalized true anomaly parameter taking into account the relevant PN contributions, yielding a 
complete description of the classical motion in time of the binary. Our main results are the explicit solution of the 
equations of motion on the one hand and the time-domain gravitational waveforms written out in detail on the other. 
To facilitate the use of our expressions, the coefficients of the waveforms are given in tabular form. We also present 
detectable waveforms for realistic sources on eccentric orbits, by numerically integrating the time dependence of the 
polarization states. As an extension of our previous work [28j we investigate the emitted gravitational wave signal 
of close hyperbolic encounters of compact objects in detail. Our main results represent the 1PN corrections to the 
expressions of Turner [23[ . These waveforms have a burstlike character and can be applied to improve present search 
algorithms for the detection of gravitational waves. 

Our paper is organized as follows. In Sec. II we describe the radial motion of the binary and introduce a suitable 
radial parametrization, the generalized true anomaly, which is valid for all type of orbital motion, namely for the ellip- 
tic, parabolic and hyperbolic cases. We analyze the connection between the coordinate time and the parametrization 
in Sec. III. In Sec. IV we describe the evolution of the polar angle, T, in the orbital plane. In Sec. V our results are 
specialized for the case of circular orbits. Section VI contains our main results, where the polarization states of the 
emitted gravitational waves are given for circular, elliptic and open orbits. In Sec. VII time-dependent gravitational 
waveforms are computed for such realistic systems as the Hulse- Taylor and the J0737-3039 pulsars, as well as binaries 
on parabolic and hyperbolic orbits, composed of 8M and 13M masses. Section VIII contains our conclusions and 
some technical details are relegated to Appendix A. Appendix B contains the constant coefficients of the waveforms. 

Throughout the paper we use units in which G = c = 1. 



II. DESCRIPTION OF THE MOTION AND THE PARAMETRIZATION OF THE ORBIT 

Similarly to [28j we introduce two coordinate systems for the description of the classical motion of the binary and the 
polarization states. The z axis of the invariant coordinate system, in which the evolution of the orbital elements and 
the dynamics of the binary are most conveniently described, is fixed to the direction of the total angular momentum 
L. The calculations can be simplified if the x and y axes of this system are chosen in a way that the direction of the 
line of sight has the components N = (sin 7, 0, cos 7) with the constant angle 7 of L and N. In the invariant system 
the separation vector r = rn is written as r = r (cos T, sin T, 0), where T is the polar angle of the motion in the 
orbital plane. 

The expressions for the polarization states become simpler with the description of the main vector quantities and 
the transverse-traceless (TT) part h^ T of the radiation field in the comoving system. In this coordinate system the 
z and x axes are aligned with the direction of the Newtonian angular momentum L^r = /ir x v and the separation 
vector r, respectively, where /1 = ra\milm is the reduced mass, m = m\ + 7712 is the total mass of the binary and v 
is the relative velocity vector. The transformation between the two systems is given by a rotation around the z axis 
with the angle T. 

The true anomaly parametrization, introduced for the description of Keplerian orbits, can be generalized to the 
perturbed Kepler problem. The form of the parametrization is determined by the radial equation of the motion, 

. 2 _ 2E 2m L 2 . 2 v 
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where rj = fi/m. The total energy E and the magnitude L of the total angular momentum are constants of motion 
which follow from a Lagrangian description [29J. The turning points, defined by f 2 = 0, are solved up to the required 
order, 
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ftfg = („ - 7)^ ± („ + 9)g T (3* - 1)^ , (2) 

where Aq — (to 2 /j 2 + 2EL 2 / /i) 1 / 2 is the magnitude of the Laplace-Runge-Lenz vector. To introduce the radial 
parametrization which is valid for all the possible values of the total energy, first we look at the following integral: 

T(r') = T(r)|- mj „= / ±f . 

J Tmirt 

The equation for the polar angle T can be evaluated following the description in [28], 

(3) 
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As a consequence of the definition T(r ro j n ) = 0. In the case of elliptic orbits r m i n and r max are the turning points 
of the radial motion, while for unbounded orbits r m i n is the only positive root, and r max — oo. According to these 
cases the value of T(r max ) is the following. For elliptic and parabolic orbits we have 

^ ij'max) | elLiptic,parabolic — 7r/C , 

where K, = (l + 3m 2 /j, 2 /L 2 ) . In the hyperbolic case 

mfiy / 2Ejl(r]L 2 E - 12/j 3 to 2 - IbLE 2 ) 



^ (j'max^l hyperbolic — Xhyp^ 
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where Xhyp = arccos (— mfx/Ao) is the maximum value of x f° r the hyperbolic motion in the Newtonian case. 

To get equivalent parametrization with the results of Damour and Deruelle [l3| and that of [29j our generalized 
true anomaly parameter x is introduced as follows 

~ir\ X— =T(r^)//C. (4) 



The solution of the radial motion is 



d(cosx) 

r(x) =r N (x) +rp N {x), (5) 



where 
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fi(fj,m + Aq cosx) 

2(6 - T])m 4 fi 6 + 2(10 - iri)EL 2 m 2 ^ + (1 - 3r])E 2 L 4 2(2 - T])mEL 2 + (6 - r/)m 3 /j 3 



(7) 



PN ^ 2A Q fi 3 (fim + Aq cosx) 2 pb{nm + A cos x) 2 ^ 

We note that the above solution is valid for all three types of orbits. 

The relation between the generalized true anomaly parameter and the coordinate time is given by 

^ = 7^ = ^r{ 1_ 2i[ (7?_ 13 ) m V + ^ - ^ A o + - 8 ) m ^° cos *] } • 

The type of orbital motion, namely elliptic, parabolic or hyperbolic, is determined by Aq in Eq.(J5]) through the 
value of the total energy E. For elliptic orbits the energy is negative and m/j > Aq. In this case the length of the 
separation vector is changing periodically between the turning points and the coordinate time t is a monotonically 
increasing function of the true anomaly parameter. For parabolic orbits the energy is zero, m/i = Aq, and r(x) 
diverges at x = 7r - bi this limit t also diverges since dt/dx ~ r 2 . For hyperbolic orbits the energy is negative and 
m/i < Aq. The divergence of r(x) as the parameter approaches the value x = cos _1 (— m/j/Ao) defines the boundary 
of the orbit. 

Because of the fact that the formal expressions of r(x) and dt/dx are identical for all types of orbits the integration 
of the angular variables formally leads to the same expressions. The difference between the three cases arises in the 
evaluation of the time dependence of x an d r, which is determined by the value of Aq. 
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III. THE COORDINATE TIME AND THE PARAMETER 



To obtain time-domain gravitational waveforms it is necessary to express the time dependence of the true anomaly 
parameter \. This time dependence, however, cannot be given analytically. The most straightforward way is the 
numerical integration of the differential equations, nevertheless, in the case of elliptic orbits the separation of different 
contributions belonging to different PN orders is highly nontrivial. As a solution of this problem first we analyze the 
analytic form of the inverse function t(x) and determine a general form of the inverse function x(t). 



A. Integration of t(\) 



In the leading Newtonian order the connection between the coordinate time and the parameter is governed by the 
relation 

dt _ jir 2 N _ L 3 



dx L fi(uifj, + A cos x) 2 
The solution of this equation for elliptic orbits is 



where 



For parabolic orbits we have 



2mfi 3 / 2 LA s'mx 
f = (-2^)3/2 Ael1 + 2E(A Q cosx + rr^) ' ( } 



^arctg (m ^ )tg ^ 
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y/m 2 /! 2 — Aq 

gtg(j) [tg 2 (f)+3] 
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while for hyperbolic orbits the result is 

_ 2m/j 3 / 2 t LA a sin x 

t — , ; „ /„ Ah 



(2S) 3 / 2 vp 2E(A cos x + mfi) 1 
where 

i hw a,) tg(f) 

\j A\ — m 2 jj, 2 

For a complete description the perturbative solution of the relation between the coordinate time and x IS required. 
Up to 1PN the corrections of r(x) m Eq.([7]) are taken into account with the neglection of higher order terms. The 
solution for elliptic orbits is 

LAosinx ma 1 / 2 rj 

1 = 2E { A Q cosxlm») ~ 4E(2E)^ ^^-IfftA* 



A cos x[4(3^ - l)EL 2 + 7(1 + ?i)m V] + m / u[4(2?j + 26)EL 2 + 7(31 + 3?y)mV] 

8A Q n 2 (m/i + A cos x) 2 

For parabolic orbits we have 

t _ [(377 - 18)m 2 M 2 + 2L 2 ] tg 2 (§) - [(9r) + 6)m 2 /i 2 - 6L 2 ] 

\2p?m 2 

and for hyperbolic orbits the solution is 



(9) 



LAasinx m/j, 1 / 2 



A cos x[4(3r] - l)EL 2 + 7(1 + ri)m V] + mfi[4(2rj + 26)EL 2 + 7(31 + 3r/)m 2 /j 

8A /j 2 (m/i + A cos x) 2 



(10) 
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B. Integration of x(f) 



The differential equation which determines the time dependence of the parameter x is 

= ±_ L , 1 

dt /ir 2 

L 2LrpN 1 



1 1 + [iV ~ 13)m V + (3ry - 1) A% + {3 V - 8)m/Mo cos X ] J 



where rjy and rpjy are given in Eq. ([5]). Since we look for the solution as x(t) = XN(t) + XpaKOi the different 
corrections to the above equation has to be separated. This process can easily be done in the case of open orbits, 
where x(t) is a bounded function at every PN order. In this case the differential equations for x are 



dXN ( dx 



dt V dt 
dxPN ( dx 



dt V dt 

' PN 

/j,(/j,m + A cosxn) 2 



_ n(pm + A cos xtv) 2 
L 3 

N 

_ 2A q xpn sin XN^rn + A cos xn) _ 2/i 2 (^m + A cos XN) 3 rpN (xn) 
L 3 L 5 

PN 
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[(n - 13)mV 2 + (3ry - 1)A 2 + (3?/ - 8)m f iA cos X n] , (12) 



where 



, , 2(6 - r])m 4 fi 6 + 2(10 - 3rj)EL 2 m 2 fi 3 + (1 - 3rj)E 2 L 4 2(2 - r])mEL 2 + (6 - ri)m 3 /j, 3 

rpN(XN) = 771 3? — A ^2 COSXJV j — j ^ ■ 

lAofi^^m + A cosxn) /.i(ptm + A cosxn) 

Since x an d xpn are bounded functions of time the equations can be integrated numerically without difficulty. 

For elliptic orbits because the time dependence of xpn is given by an unbounded function, therefore simple series 
expansion gives rise to secular divergences. To deal with this problem one has to separate the bounded and unbounded 
terms in xn and xpn- 

To find the solution we note that t(x) can be written as 

t(x) - Tx + f{x) 

= t N ( X ) + t PN ( X ) - (T W + T PN )x + f N (x) + fPN(x) , (13) 

where Tn and Tpjsr are constants, and the functions Jn(x) and /pjv(x) are bounded. To prepare the numerical 
inversion of i(x)> first we use the following general theorem about smooth functions. 

Let F(x) be a smooth function, which can be written in the form y — F(x) = Ax + B(x), where A ^ is constant, 
and B(x) is a bounded function of x. If the F^ 1 inverse function exists, writing it in the form F^ 1 {y) = y/A + b(y), 
b(y) will be a bounded function of y. 

With the use of this theorem we look for the solution in the following form: 

x{t) = Qt + g(t) 

= n N t + g N {t)+n PN t + g PN (t) . (14) 

The constants TV, Tp n, and hence Qn and flpN can be determined by the investigation of the circular orbit limit. 
The results are 

{-2E) 3 ' 2 ( V -15)E 3 

— T77T~ , "PAT 



m/i 3 / 2 ' ^" (-2^) 1 /2 TOAt 5/2 ■ 

Moreover, we note that the use of the straightforward series expansion during the calculation of the sine and cosine 
of x(t) results in the appearance of secular divergences. To solve this problem, as an example, the following expression 
is to be used: 

sinx(f) = sin [fijyf + flpNt + ffjv(t)] + cos [fijvf + QpNt + gN(t)]gpN(t) , 

and similarly for cosx(t), which leads to a combined amplitude and frequency series expansion. In our calculations 
we use the above expansion to avoid secular divergences in the description of the motion of the binary and the time 
dependence of the polarization states of the detectable waveform. 

Although the above expansion helps us to avoid the appearance of secular terms, it gives rise to technical difficulties 
when separating the different PN corrections on the level of Eq. pip . With the integration of both the Newtonian 
order and the full equation, one can determine <?jv(t) and gpN(t) separately, and give the full description of the time 
dependence of x- The details of the calculations are given in Appendix A. 
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IV. ANGULAR EVOLUTION 



Similarly to r the polar angle T is decomposed as 



T = T N + T PN . 

To solve the equations of motion the form of the relative velocity vector is required in terms of the Euler angles. In 

the invariant system it is expressed as v = fr, rT,0j. For the description of gravitational waves and the equations 

of motion one needs an alternative form of the components of the velocity vector, both in their general form and in 
terms of the parametrization. 

In terms of the constants of motion one can express f 2 , Eq.([I]), and the square of the velocity vector, which is 

fF,\ 2 

3 



2E 2m 



Kl-*>)(f) + 2(6 -7,)^ + (10- 5,) (^) : 



ijmL 2 



H r 

From our previous results if follows that v?, — r 2 and, moreover, 



2 2 -2 

v = v — r 



L 2 , . EL 2 , ,mh 2 
( 2 - 6 ^-( 8 - 4 ^ 



H 2 r 2 



This indicates that 



L 



L 



vj_ = 



r PN „ s E , , n . m 

— + (1 - 3ry)- + (4 - 2ry) — 
r N [i r N 



When we insert the parameter dependence of r/v and rpN in the above formulas the components of the velocity vector 
become 



and 



A)Sinx 2(77 - l)m 4 ^ 6 + 3(3?7 - 1)£; 2 L 4 + 2(4r/ - h)EL 2 m 2 ^ . (8 - 3ry)m^ 

fim + A cosx 2(9 - 1077)771 V + 2(7 - ?,rj)EL 2 ^m 2 + 3(1 - 3ri)E 2 L 4 



L 



2fi 2 A Q L 3 



cosx + 



(2 + jj)m 3 fi 3 + (3 + r))mEL 2 2(r\ - 2)m 3 /j 3 + 4(77 - 2)mEL 2 



2L 3 



From the relation tii = rT we have 



T = 



1 - 



I? 



2rpN , , E 777 

+ (1 - 37?)- + (4 - 27?) — 
r n A* Tn. 



cos x , 



(15) 

(16) 
(17) 

(18) 



which can be expressed in terms of the parametrization as 



/ dT \ 6m 2 /_i 2 + r]mfiA cos \ 

\dx~J = 2V^ 



The solution for T is 



T = T + x 



6m ju x + 77777/iAo sinx 
2L2 



(19) 



The separation vector and the waveforms contain expressions of sinT and cosT, and their higher harmonics. With 
the use of a series expansion, for example, sin T can be expressed as 

6m 2 fi 2 x + 77777 {j,Aq sinx s 



sin T = sin (T +X) ~ cos (To + x) 



2L 2 



For elliptic orbits there is no upper limit for x as a function of the coordinate time and the (3m 2 /a 2 /L 2 ) xcos (To + x) 
perturbative term gives rise to a secular divergence. To handle this problem we follow the method described above 
in the case of the time dependence of the parameter, namely, instead of the above expression we apply the following 
decomposition: 



sin T = sin ( To + x 



3m n x 



cos T + x 



3m fj, x \ T] m ^AQ smx 



L 2 J \ u A L 2 J 

where no linear terms appear (note that a similar expression is valid for cosT). 



2L 2 
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V. EQUATIONS OF MOTION FOR CIRCULAR ORBITS 

The relative velocity vector can generally be written as v = fn + ram, where m is the unit vector in the direction 
of the y axis of the comoving system. The definition of a circular orbit is 

r = uj = . 

From the above decomposition of the relative velocity it turns out that u) = T. As a consequence, 



v = u± 



L 
///• 



, E , . m 

1- 1-37? -- (4-277 - 
fi r 



(20) 



In the circular orbit limit the constant radius r of the orbit cannot be evaluated analytically in all orders of the 
approximation and hence in our calculations we keep r as a constant parameter determined later. 
From Eq.(|20p the time evolution of T can be written as 



T 



Lt 
fir 2 



1 - (1-37/)- - (4-2??)- 
fi r 



(21) 



During the evaluation of quantities including the sine and cosine of T there are difficulties occurring in connection 
with the secular divergences as in the elliptic case. The solution of this problem is the same, although it leads to 
simpler form. For example, 



sin T = sin < — ^ 
I A" 



1- (1-37?)- -(4 -27?)- 
/i r 



In this case the secular terms can be treated with a simple expansion series of the orbital frequency. 

VI. THE STRUCTURE OF THE DETECTABLE WAVEFORM 

The parameter and, for circular orbits, time dependence of the waveform and the polarization states h+ and h x are 
determined by the use of the method described in [28J . The polarization states of the detectable gravitational waves 
are calculated by the projections of the transverse-traceless tensor h l ^ T representing metric perturbations, 



■^(PiQj + UPj)h 



TT ) 



where p is a unit vector in the orbital plane perpendicular to the direction of the line of sight N, and q = N x p. 

Since we have chosen the comoving system to describe this projection of hj, T we determine the components of N, p 
and q in this system. With the use of the transformation law between the two coordinate systems N has the following 
form in the comoving one: 

sin 7 cos T 
N = | — sin 7 sin T 

cos 7 

Moreover, the conditions for p determine its components as 




To separate the different contributions the transverse-traceless part of the radiation field h l ^ T can be decomposed in 
the post-Newtonian approximation as [30l | 
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where D is the distance between the observer and the source and we have collected all the terms which are relevant 
up to 1 PN order. Q IJ denotes the quadrupole (or Newtonian) term, P 5 Q^, and PQ y are higher order relativistic 
corrections. The detailed expressions for these contributions are given in 130, 311. 
We choose a similar decomposition for the polarization states h+ and h x , 

h N+ +h 0.5+ +h l+^ ( 

It is worth mentioning that this decomposition is formal in the sense that higher order terms in the description of 
motion can cause higher order contributions in h * . 

Since the description of the motion of the binary system has different properties for the different types of the orbit, 
and hence one can use different simplifications throughout the evaluations, the elliptic, circular, and open orbit cases 
are discussed separately including the description of the simplifications and the structure of the detectable waveform. 
The constant coefficients are given in detail in Appendix B. 




A. The elliptic orbit case 

For elliptic orbits a combined series expansion has to be included to deal with the problem of secular divergences 
described above. These secular terms give rise to a frequencylike series beside the usual amplitude expansion. In 
order to define this we introduce 

, 3mV , 00 x 
X = T + X J2~ x > ( 22 > 

which is valid up to 1PN. To lowest, Newtonian order the above definition becomes simply \' — To + X- Since our 
goal is to give waveform expressions valid up to 1PN in the following in our expressions the original definition of x' 
will be used. 

Following the main steps described above, the parameter dependence of the detectable gravitational waveform for 
eccentric orbits can be given as 

^=EE [c"}co S (k X ' + jx) + s"* S in(kx'+jx)] , 

fe=-2 j=-3 



h ° K = E c °s + JX) + < ■ x sin (k X ' + jx)] , 

k,j=-3 



hlt = E [cl t j ™s(kx'+3X) + s£ j sm(k X '+jx)] , (23) 

k,j=-4 

where the S and C coefficients are depending on the energy E, the magnitude of the total angular momentum L, the 
masses of the objects and the 7 angle, and hence they are constants. The explicit form of these constant coefficients 
is given in Appendix B. 

It is worth mentioning that the above expansions are formal, since the time dependence of the true anomaly 
parameter \ contains higher order corrections. These can modify the expressions for both the secular terms and the 
amplitude expansion. Completing the analysis one arrives to a combined frequency and amplitude series expansion. 
In this work the whole method is used to evaluate (numerically) the time dependence of the polarization states shown 
in the next chapter. 



B. The circular orbit case 



As it has been shown in the circular orbit case we use the r constant length of the separation vector as a parameter 
throughout our calculations. In this case Eq. (fiZ!?j) reduces to 

X' = To + — t - -A [(1 - 3v)Er + (4 - 27 1 )m^} = T + u N t + u PN t . 
fir ji v 
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With the use of this notation the structure of the detectable waveform becomes 

hNt = E ^ cos + < X+ sin (fcx')] , 

fe=0 



^ = E [^°' 5x cos (fcx') + ^°- 5x sin (k X ')] , 

k=0 



h 1 * = E cos + sin (**')] • (24) 

fc=0 

Again, the S 1 and C constant coefficients are depending on E, L, the masses, r, and the angle 7. The explicit form of 
these coefficients is given in Appendix B. 

The high frequency terms in h N arise from the 1PN terms of the description of motion. At the Newtonian order the 
frequency of the detectable waveform is twice the orbital frequency ujn, in agreement with the quadrupole formalism. 

C. Waveforms for open orbits 

The evaluation process of the waveform is almost the same as it is in the elliptic case. The only difference is that 
there are no secular terms arising in the parabolic and hyperbolic cases, and it is possible to use the usual expansion 
of the sine and cosine of the angular variables. The parameter dependence of the polarization states can be written 

as 

3 + + 4 + + 

^=XE [°kx cos + S kx sin ( fc x)] + E i C k X cos ( fc x) + sin (kx)] , 

fc=0 fe=0 



h° K = E [°1' 5X cos + S fe °- 5x sin {k X )] , 



fc=0 

where the constant coefficients are depending on the same parameters as in the elliptic case, and, moreover, on the 
T initial value of T. We note that the first type of terms in h N are linear in the parameter. In these terms all the 
coefficients are at IPN order, and since the value of the parameter is bounded, the description is self-consistent. 

VII. REALISTIC WAVEFORM SIGNALS 

Following the description above, after the integration of the time dependence of the generalized true anomaly 
parameter the realistic, time-dependent gravitational waveforms can be given. In this section we present Newtonian 
waveforms as well as 0.5PN and IPN corrections for well-known or reference systems in order to analyze the effects 
of the orbital eccentricity. 

A. Waveform of realistic elliptic binaries 



For further investigation, we numerically integrate the time dependence of the polarization states of the detectable 
waveform for two realistic and well-known sources, namely, the Hulse- Taylor and J0737-3039 pulsars, see Table I. Our 
goal is to investigate the effects of the eccentricity of the orbit and the 0.5PN and IPN corrections to the waveform. 
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Parameters 


Hulse- Taylor 


J0737-3039 


Mass of object 1 (mi) 


1.441 M 


1.337 M 


Mass of object 2 (ma) 


1.387 M 


1.25 M 


Eccentricity (e) 


0.617 


0.0878 


Periastron (a) 


1.950M0 9 [m] 


0.8787- 10 9 [m] 


Distance from detector (D) 


21000 [ly] 


1956 [ly] 


Orbital frequency (Torb) 


27907 [s] 


8834 [s] 


J-N angle (7) 


21.12° 


28.89° 



Table I: Parameters of the investigated sources 




Figure 1: Time dependence of the polarization states (h+ depicted by solid line and h x corresponds to dashed line) at the 
Newtonian order for the Hulse- Taylor (left) and for the J0737-3039 pulsars (right). The amplitudes are multiplied by the factor 
10 20 . 

To keep our figures transparent, the h+ and h x polarization states will be shown on the same plot. We collect 
the results belonging to different sources and PN orders. This simplification allows us to show the main differences 
between the selected systems order by order. 

Following these guidelines, first we show the polarization states at the Newtonian order. 

Figure [1] shows that the eccentricity of the orbit leads to many differences between the waveform predictions of the 
different sources. In the case of the Hulse- Taylor pulsar, the results for the emitted waveform are significantly different 
from the usual circular orbit case, since the eccentricity gives rise to many higher harmonics of the waveform. In the 
case of J0737-3039 pulsar, the eccentricity is almost 1 order of magnitude smaller and the emitted waveform is similar 
to the one generated by binaries on the circular orbit. We note that these differences are much more significant than 
the higher corrections of the PN series in both cases. 

The amplitude of the emitted waveform of J0737-3039 is much higher than that of the Hulse- Taylor pulsar since 
the periastron of the former is less than half of the latter one, and it is much closer to the detector. 

The higher order corrections show the same properties of the waveforms for the different sources. It is worth 
mentioning that in the circular orbit limit the frequency of the emitted waveform is twice the orbital frequency of the 
motion. 




Figure 2: Time dependence of the first corrections to the polarization states at 0.5PN order. Again, h+ is depicted by the solid 
line and h x by the dashed line. The corrections for the Hulse- Taylor pulsar are shown on the left and for the J-0737 pulsar on 
the right. For further comparison amplitudes are multiplied by the factor 10 20 . 
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Figure 3: Multipolar corrections to the polarization states at 1PN order. Again, the corrections for the Hulse- Taylor pulsar 
are shown on the left, and for the J0737-3039 pulsar on the right. Solid line shows the h+ polarization state, and dashed line 
shows h x - Again, the amplitudes are multiplied by the factor 10 20 . 



3. x io -7 r 




Figure 4: 1PN order corrections to the polarization states, which arise from the perturbative description of the elements of the 
motion. The corrections for the Hulse- Taylor pulsar are shown on the left and for the J0737-3039 pulsar on the right. The 
solid line shows the h+ polarization state and the dashed line shows h x - To compare the different terms of these corrections, 
the amplitudes are multiplied by the same factor as before, namely 10 20 . 




Figure 5: In these figures we plot the full 1PN corrections to the polarization states for the Hulse- Taylor pulsar (left) and for 
the J0737-3039 pulsar (right). The solid line represents the h+ polarization state and the dashed line shows h x - 

Beyond the properties described above there are two additional effects on the waveform shown by Figs. [2KSJ By 
comparing the leading amplitudes of the different corrections it turns out that, since the masses of the objects are 
almost equal for both sources, and the 0.5PN order corrections are multiplied by Sm, these corrections to the waveform 
are suppressed by the simple 1PN corrections. 

Since the lowest order corrections to the dynamical variables arise at 1PN, the 0.5PN corrections to the waveform 
are originated only from the multipolar expansion. For 1PN corrections both the multipolar expansion and the 
contributions from the description of the motion contribute to the waveform. On the other hand 1PN corrections 
from the orbital motion of the binaries are much more significant than the contributions from the multipolar expansion 
of the post-Minkowskian series. 

Although a full analysis of the waveform would require Fourier expansion, the main properties of these realistic 
waveforms are transparent. The contributions which can be seen only on long-term calculations are from the 1PN 



12 



corrections to the frequency of the waveform arising through the time dependence of true anomaly parameter and the 
periastron precession. These effects cause the slow shift of the basic frequency of the waveform. 

B. Waveform of reference parabolic and hyperbolic orbits 

Our main goal is to evaluate the first corrections to the detectable gravitational waves emitted by binaries scattering 
on hyperbolic and parabolic orbits. 

To investigate the properties of the emitted waveforms of scattering binary systems, we introduce two hypothetical 
reference binaries, see Table II. The masses of the orbiting bodies are chosen following the parameter estimations in 
0], although we determined the values in a way which avoids the suppression of the 0.5PN corrections. The minimal 
distance between the objects and the distance from the detector are chosen to be the same as in the case of the 
Hulse- Taylor pulsar. It will give the opportunity to compare the waveform properties between open and closed orbits. 

To investigate sources in which the emission of GWs results in the formation of a closed binary (in extreme cases 
even zoom-whirl orbits) we discuss parabolic sources. To analyze the effects of the eccentricity, besides parabolic 
orbits we introduce a hyperbolic source (e = 2). 



Parameters 


Hypothetic sources 


Mass of object 1 (mi) 
Mass of object 2 (rn.2) 
Eccentricity (e) 
Minimal distance (r min ) 
Distance from detector (D) 
J-N angle (7) 


8 M 
13 M Q 
1 and 2 
1.950M0 9 [m] 
21000 [ly] 
45° 



Table II: Parameters of the investigated open sources 




-0.1 L -0.1 - 



Figure 6: Time dependence of the polarization states at Newtonian order for the two open orbit sources. The solid line shows 
the parabolic case and the dashed line shows the waveform of the hyperbolic (e = 2) source. The h+ polarization state are 
depicted on the left figure and h x is shown on the left one. The amplitudes are multiplied by the factor 10 20 . 

The waveform for open orbits is not periodic (or quasiperiodic), it is burstlike, as it is shown in [23l |. Our hypothetic 
sources provide a high amplitude signal in the scale of a few hours. Binaries in Q provide gravitational waves with 
higher amplitude, since the minimal distance between the objects are much smaller than ours. 

In Figs. IBTfTOl it can be seen that the eccentricity of the orbit affects the frequency of the emitted signal. With the 
same properties, the higher the eccentricity, the more spiky the signal with higher amplitude at each PN order. 

VIII. CONCLUDING REMARKS 

We have introduced a generalized true anomaly parametrization for binary systems on either closed or on open 
orbit and applied it to evaluate gravitational waveforms emitted by them. In addition to elliptic orbits we consider 
the scattering of binary systems following quasiparabolic and quasihyperbolic orbits. Our new parametrization leads 
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Figure 7: Time dependence of the first corrections to the waveform polarization states at 0.5PN order. Again, the corrections 
to the h+ polarization are shown on the left and for h x on the right. The solid line shows the parabolic case and the dashed 
line shows e = 2 hyperbolic source. For further comparison amplitudes are multiplied by the factor 10 20 . 
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Figure 8: Multipolar corrections to the polarization states at 1PN order. Again, the corrections to h+ are shown on the left 
and to h x on the right. Parabolic and hyperbolic orbit cases are depicted by solid and dashed lines, respectively. Again, the 
amplitudes are multiplied by the factor 10 20 . 
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Figure 9: 1PN order corrections to the polarization states, which arise from the perturbative description of the elements of 
the motion. The notation is the same as on the former figures in this subsection. To compare the different terms of these 
corrections, the amplitudes are multiplied by the same factor as before, namely 10 20 . 



to a universal description valid for elliptic, parabolic and hyperbolic cases. The only differences are in the numerical 
integration determining the time dependence. With the use of this parametrization the radial and angular equations 
of the motion are solved, and with the inclusion of the results the analytic form of the polarization states for the 
emitted gravitational wave signals are evaluated for all types of sources including circular orbits, too. 

The time dependence of the generalized true anomaly parameter x is determined numerically for arbitrary parameter 
values. We have evaluated polarization states of the emitted radiation for realistic, and some well-known astrophysical 
sources. For elliptic orbits we presented the corrections to the waveform order by order for the Hulse- Taylor and J0737- 
3039 pulsars. For open orbits we have computed the waveforms for two sources with physically reasonable orbital 
elements to illustrate their main characteristics. 
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Figure 10: In these figures we plot the full 1PN corrections to the polarization states for the parabolic (soled line) and hyperbolic 
(dashed line) sources. The corrections to the h+ polarization states are depicted on the left, while to h x on the right. 



Our results illustrate the importance of the effects of eccentricity arising from the presence of higher harmonics in 
the waveform even at Newtonian order. The eccentricity increases the possibility of the detection of GWs emitted by 
sources on such orbits even if the eccentricity is very small, as it is in the case of the pulsar J0737-3039. 

Our investigation of open orbit sources show that they provide unique, burstlike wave signals, and this behavior 
becomes more characteristic for higher order corrections. Meanwhile, these sources and the emitted waveforms can 
be described analytically, only the time dependence of the parameter has to be treated numerically. It gives the 
opportunity to investigate and create data banks for these sources to use matched filtering techniques. 

Though the event rate estimations show that the detection of hyperbolic orbits are less probable, the emitted signal 
is more "spiky" for hyperbolic orbits than for parabolic ones with higher amplitude and frequency. These differences 
become more significant at higher PN order. We have chosen a quite large impact parameter (namely, almost equal 
to the periastron of the Hulse- Taylor pulsar) to compare the differences between closed and open orbit cases, and to 
investigatejhe effects of the eccentricity on the widest possible range. With the application of the parameter values 
chosen in Q, the frequency and the amplitude of the emitted waveform significantly increases. 

The burstlike behavior of the waveform for open orbits, and the possibility of an analytic description, make these 
systems an important source of detectable gravitational waves for LISA. The significantly smaller impact parameter 
choice makes these systems potentially detectable for even the advanced ground-based detectors. The higher order 
description gives rise to many additional advantages, since it provides higher frequency terms for the matched filtering 
evaluations without the assumption of extreme parameter values for the astrophysical sources. 

The results detailed above provide the possibility to generalize this formalism to describe spin and quadrupole 
effects on the detectable waveform. Furthermore, at higher orders it is possible to describe the closure of the orbit 
controlled by radiation reaction. 
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APPENDIX A: NEGLECTING INTERFERENCE TERMS IN X (i) 

When solving Eq. lfTTj) for elliptic orbits, the Newtonian and perturbative terms cannot be straightforwardly sepa- 
rated because of the secular divergences. The easiest way to solve the differential equation numerically is to solve the 
Newtonian and 1PN order equations with either the integration of the differential equations or with the inversion of 
the analytic solution t(x)- 

When the different order contributions in the resulting numerical data are to be separated it is nontrivial to neglect 
the interference terms. In the following we give the description of the problem, and give the guidelines for the solution. 
The Newtonian order equation has the following structure 

-T7- = Mcosxat) , (26) 
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and similarly the 1PN order equation is 



dx 
dt 



= /(cosx) , 



(27) 



where 

/(cosx) = /w(cosx) + /p(cosx) , 
and we wish to find the solution of the form 

\N = &Nt + QN(t) , X = i^N + Qp)t + g(t) 



(28) 



Using the above results one can evaluate both the gN{t) and g(t) functions (SI n and Sip are known from analytic 
evaluations, see Sec. III.B.). To be able to determine the contributions of the different order to the 



X = (fijv + ttp)t + g N (t) + g P (t) 

solution, our main goal is to determine the function g]\r(t). 
To achieve this we rewrite the above differential equations as 



dgN(t) 
dt 



/at (cos [Sl N t + g N (t)}) - 



and 



d[g N (t) + g P (t)} 
dt 



= f N (cos [(Ojv + Q P )t + g N (t)} - g P (t) sin [(Sl N + Sl P )t + g N (t)}) 
+ / P (cos [(fijv + Sl P )t + g N (t)} - [Sl N + Sl P ) + o(e 2 ) . 



After introducing the 



, S]w — Sip 
t' = —t 



variable we will analyze the equation 
d{g N {t') + gp{t')] 



dt 



This equation will become 



= f N {cos[fl N t + g N (t')]-g P (t')sm[n N t + g N (t')]) [ 1 
+ / P (cos [Sl N t + g N (t')] -Sl N + o(e 2 ) . 



si P 



d~g N {t') 



dt 



= f N (cos[Sl N t + g N (t')]-Sl N 



N 



at the Newtonian order. Comparing it with Eq. (|30p it turns out that 



git') = 3Jv(i) = .9« 



1 



n 



N 



With the use of it the whole solution of the problem is 

Sl P 



X = (Ojv + £lp)t + g N 



1 



si 



N 



t) + g(t) - g N 



si~ 



N 



where gN(t) and g(t) are to be determined numerically, SIn and Sip are evaluated analytically. 



(29) 



(30) 



(31) 



(32) 
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APPENDIX B: CONSTANT COEFFICIENTS FOR THE WAVEFORM EXPRESSIONS 

Our analytic calculations give rise to complicated expressions for the polarization states of the gravitational wave 
signals in all relevant cases. In the previous sections the main structure of the waveforms was given and their main 
properties were analyzed. In this Appendix we collect the explicit form of the constant coefficients appearing in the 
waveform expressions. 

We start by introducing the following short-hand notations: 

mil 

K = 1 — 3?7 , B = 8 - 3ry , D = — - , 

Aq 

A N _ A A N _ A m ^ A N _ TO V 

A i ~ 4£a ' A 2 - -§I2- > A 3 - "^r , 

N _ A 3 mfj, JV _A|mV a n _ ApmV w _ _E^_ _ _ mV 

4 ~ 8L 4 ' 5 ~ 16L 4 ' 6 ~ 4L 4 ' ^ 7 " 16 M 2 ' ^ 8 " 8L 2 ' 9 " 16L 4 ' 

A N - DA N A N - HA N A N — DA N 

^10 — U-<^7 > -fill — -^^S ' -™-12 — i 

4 0.5 ^mAg 4 o.5 _ SmAlu 5 _ <5mA TO/i 2 5 _ 5mm 2 /i 3 
1 " 256L 3 m ' 2 ~ 128L 3 ' 3 ~ 64L 3 ' 4 ~ 32L 3 ' 

Al = Sa ' ^ 2 = ^ /384 ' A " = A " /48 ' Al = A " m ' ^ = A * /3 ■ (33) 

Although the above coefficients are not independent, they are introduced separately to keep our tables clear and 
well-defined. 

As it can be seen from the results of Sec. VI. in the analytic waveform expressions, Eqs. (l23|) . the coefficients Cf 
and Sp of \ an d x' are time independent. These coefficients depend on the different harmonics of the angle 7 and 
we introduce the following general decomposition valid at least up to 1PN order and for both polarization states: 

00 

C f = Yj [ a ? F sin (h) + b ?" cos (h) , 
1=0 
00 

S F = Y^ [ a f F sin ih) + b i F cos (h) ■ (34) 
i=o 

Here F is a multi-index containing the index of the polarization state, the PN order of the expression, and the index 
of the harmonics of x and/or The coefficients af F , bf F , af F ,and bf F in Eq. (f3"4"|) will be given in tabular form 
for both polarization states and PN orders. Any coefficients af F 7 bf F 7 af F and bf F not listed in the tables are 
identically zero. With the use of these coefficients listed in the tables the quantities Cf and S F in Eq. (|34[) are easily 
reconstructed. 

A. Circular orbit case 

Let us recall the results presented in Sec. VLB. The polarization states for the detectable gravitational wave signal 
can be written as 

2 + 4 



h N i = J2 [Ck cos (fcx') + S fc x sin (k X ')] , (35) 

fc=0 

h 0.5t = £ [ C ° k K cos (k X >) + S° k 5t sin (k X ')] , (36) 

k=0 



k=0 



where 



X' = T + - -^3 [C 1 - ^)Er + (4 - 2t?)to M ] 
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Instead of using the short-hand notations of Eq. 

L 2 



a?' 5 



4r 2 /i 2 ' 
SmL 3 
32mn 3 r 3 

L 4 
38V 4 r 4 ' 



m 
4r ' 



, we use the following expressions: 
EL 2 



3 =(1-3t?) 



N 



2r 2 /x 3 



L 2 m 
r 3 fj, 2 



a° 2 5 



SmL 
~ 32/jr 2 

L 2 m 
384/j 2 r 3 



384r 2 



As an example, the Newtonian expression for C 2 arising with the use of Table Hill as 



a 



N+ 



cos 2j + 1>, 



N 



cos27 + 3(a 1 +a 



The relevant coefficients for circular orbits are listed in Tables Mil - IVIII 



State 


+ 


X 


Harmonic 


cos 27 


1 


cos 7 




aF + aZ 


3(af + a 2 v ) 





riN 




a? -eg 





Q N 

<J2 








4(^ + 0?) 



(38) 



Table III: This table contains the relevant bf k and bf k coefficients for evaluating every C k and S£ in Eq. (|35[) restricted to the 
Newtonian order. Any other coefficients are zero. The column under the sign " + " contains the coefficients for h N+ , and under 
the sign "x" the coefficients for h Nx can be found. 



State 


+ 


X 


Harmonic 


sin 37 


sin 7 


sin 27 




-4a?' 5 


-20a?' 5 


8a?' 5 - 28a?- 5 


C 0.5 


4a?' 5 


20a?' 5 


-8a?' 5 - 20a ( 2 K5 


0O.5 
*3 


-2a?' 5 + 7a° 2 5 


5(-2a?- 5 + 7aip) 


-16a° 5 


gO.S 


n 0.5 1 n 0.5 

oa! + 7a 2 


-2a?' 5 + 19a° 2 - 5 


16a? 5 



Table IV: This table details the relevant coefficients for evaluating every C k and in Eq.(j36} for both polarization states. 



The coefficients not included are zero. Again, the column under the sign " + " contains the coefficients for h 
sign "x" the coefficients for h°' 5x can be found. 



and under the 



Harmonic 


cos 47 


cos 27 


1 


cl 


K{6a\ + 51al + 7al) 


4K(6a\ + 51al + 7a\) 


5K(6a\ + 51al + 7a\) 


cl 


4K(-6a\ + 3al + 7a\) 


16[3Ka\ + 9(1 + K)a\ - 29a\\ 


4[42Ka\ + 3(7 + 87r))al + (-355 + 2177)03] 


cl 


3K(6a\ - 13a\ + 7al) 


V2[-QKa\ + (23 - 17 V )al + (-41 + 7t])ah] 


3[18Ka{ + (-79 + 29r7)a^ + (157 - 777)03] 



Table V: This table details the relevant coefficients for evaluating every nonzero Cl in Eq. (|37|l for h 1+ . It is worth mentioning 
that in this case all S k are zero. 



B. Elliptic orbit case 

Here we recall the results discussed in Sec. VI. A to collect all the expressions for the elliptic orbit case: 

2 3 

^ = EE [C^cos(kx'+jx) + S^sm(k X '+jx)} , (39) 
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Harmonic 


cos 37 


cos 7 


si 


4K(6a\ + 51a^ + 7a\) 


4K(6al + 51a^ + 7a\) 


sl 


&K(-&a\ + l5al + 7a\) 


8[30Ka\ + 3(11 + 39?7)aJ + (-239 + 2177)03] 



Table VI: This table details the relevant coefficients for evaluating every nonzero S\ in Eq. (|37p for h lx . It is worth mentioning 
that in this case all C\ are zero. 



State 


+ 


X 


Harmonic 


cos 27 


1 


cos 7 


cf 


af + af - (af + af ) 


3(af + af ) - 3(af + af ) 





r<N 
<^0 


-a? + a? + (a£ + a?) 


af - af - (af + af ) 





qN 
<->2 








4(af + af ) - 4(af + af ) 



Table VII: In this table the full coefficients are presented for Eq. (|35[) . including 1 PN order terms from the corrections of the 
description of the motion. The table is organized in the same way as Table Hill 



h°- 5 t = cos + 3X) + ST/ sin (k X ' + iX)) , (40) 

k,j=-3 



where 



h 1 * = J2 cos ^ + 3X) + Sl* sin (k X ' + 3X)] , 

k,j—-4 



1 ™ 3m 2 /1 2 
X = T + X jp— X ■ 



(41) 



Again, to collect the explicit expressions for the Ckj and Skj constant quantities we will use the short-hand 
notations and definitions of Eq. (|33[) and Eq. (|34[) . 

For elliptic orbits our tables are organized in the same way as in the circular orbit case. As an example, we detail 
the expression of C\\ with the use of Table [Xj 



Cl+ = 64^ COS [4 7 ] + % iA cos [2 7 ] 

= {K81A\ + KAh2A\) cos [4 7 ] + {i2AKA\ + 1808KAI) cos [2 7 ] - (A0hKA\ + 2260^ A\) 
The relevant coefficients for elliptic orbits are listed in Tables IVIIII - IXIIII 



r i+ 



r 1+ 



(42) 



State 


+ 


Harmonic 


cos 27 


1 


<^2,1 


Af 


3Af 


r^N 
( -'2,0 


Af 


3Af 


<^2,-l 


5Af 


15Af 


( -'2,~2 


Af 


3Af 




-Af 


Af 


r^N 
^0,0 


-Af 


Af 



State 


X 


Harmonic 


cos 7 


C JV 
"->2.1 


4Af 


C JV 
'-'2,0 


4Af 


C JV 
"->-2,2 


-4Af 


C JV 
"5-2.1 


-20Af 



Table VIII: This table contains the relevant b t h ' 3 and b t k ' 3 coefficients for evaluating every non-zero Cgj and Sgj in Eq.j39]) 
restricted to the Newtonian order. The left table contains the coefficients for h N+ , and in the right one the coefficients for h Nx 
can be found. It is worth to mention that for the "plus" state every S*f , and for the "cross" state every Cf are zero. 



Tables IXIll and IXIIII contain the 1PN corrections to the Newtonian coefficients given in Table ELLA from which the 
full coefficients in Eq. (|39[) can be reconstructed. 
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State 


+ 


X 


harmonic 


sin 37 


sin 7 


sin 27 


^.0.5 

°3,3 


-4A?- 5 


-20A?- 5 


-4A°- 5 


W).S 
°3,2 


-12A§- S 


-60A^ 5 


-4Ag 5 


/~t0.5 
°3,1 


-4(3A?' 5 + A3 6 ) 


-4(15A?- 5 + A° 3 5 ) 


4(11A?- 5 -2Ag' 5 ) 


W).5 
°3,0 


-4(6A% 5 + A * 5 ) 


-20{6A° 2 5 + A°/°) 


20(2^' 5 - A° 4 - 5 ) 


rio.5 

°3,-l 


-12(A? 5 + A° 3 5 ) 


-60(A ( l 5 + A% 6 ) 


4(13A?' 5 - 1AAY°) 


r<o.5 

°3,-2 


-12Af s 


-60A«' 5 


-60A^ 5 


^3,-3 


-4A?- 5 


-20A?' 5 


-28A?' 5 


°1,3 


4a;- 5 


4A?' S 


-4a;- 5 


^0.5 
^1,2 




6(L4§- 5 


-204 f 2 '' 5 


W).5 
^1,1 


UiAl' 5 + A° 3 - 5 ) 


60(A?' 5 + A% 5 ) 


-20(A?' 5 + 2^' 5 ) 


°1,0 


4{6A° 2 5 + A /') 


20(6A 2 1 - 5 + ^- 5 ) 


-4(18^- 5 + 


,^0.5 
^1,-1 


12(A?' 5 + A 3 ] - 5 ) 


60(A? 5 + AV°) 


-4(3A?' 5 + 22A 3 ] - 5 ) 


r<o.5 

^1,-2 


12A° 2 5 


60A° 2 5 


-84A% B 


,^0.5 


AA'l 5 


2(L4?' 5 


-28A?' 5 


c 0.5 
>->3,3 


A 0.5 




-WAf 5 


o0.5 
"3,2 


-A° 2 5 


-5A° 2 5 


-48A% B 


c 0.5 
"53,1 


-11< 5 + 2A° 3 5 


-5(11A?' 5 -2^' 5 ) 


-48(A?' 5 + A° 3 5 ) 


o0.5 
"3,0 


-5(2A° 2 5 - A / ) 


-25(2A§- 5 - 


-W(6A° 2 5 + A° 4 - 5 ) 


o0.5 
■Jl,3 


llAf 5 


-25A?' 5 


16A1 5 


o0.5 
"1,2 


23A^ 5 


-29A§' 5 


48A% S 


q0.5 
"1,1 


-23A?' 5 - 22Ag' 5 


-29A?' 5 + 14,4§- 5 


48(A?' 5 + A9/°) 


o0.5 
"1,0 


38A§' 5 + l3A ( i 5 


-3QA% 5 - 17A° 4 - 5 


16(6A 2 >5 +A° 4 - 5 ) 


aO.5 
"3-1,3 


-13A?' 5 


-17A?' 5 


-16< 5 


o0.5 
"-1,2 


~39A° 2 5 


-51A|]- 5 


-48A2' 5 


oO.S 
"-1,1 


-A ( l 5 - 93/2A% 5 


-21A1 5 - 50A% 5 


-16(A?' 5 + Ag' 5 ) 


o0.5 
"-3,3 




-35A?' 5 


16A?' 5 


c0.5 
"-3,2 






48vl^ 5 


o0.5 
"-3,1 


13A?' 5 - 14A§ 5 


5(13A?' 5 - 14A§' 5 ) 


48(A?' 5 + Af s ) 



Table IX: This table contains the relevant a l k ' 3 and a l k ' 3 coefficients for evaluating the explicit form of Eq. (l40[) . Any other 
coefficients are zero in this equation. The columns under the sign " + " contains the coefficients for h°' 5+ and under the sign 
"x" the coefficients for /i°' 5x can be found. 

C. Open orbit case 

For open orbits the structure of the polarization states of the gravitational waves can be written as 

3 + + 4 + + 

hNt = X E [°kx cos (**) + S kx sin ( fc *)] + E [ C fc * cos ^) + S k X sin (**)] , (43) 

fc=0 fc=0 

feO-Bj = J2 [C° k 5+X cos (k X ) + St 5t sin (k X )] , (44) 

fc=0 

7 + + 

k=0 

In this case there are no secular terms in the approximation scheme, and hence a more straightforward series 
expansion can be given. It makes the structure of the waveform expressions simpler, but on the other hand the 
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Harmonic 


cos [47] 


cos [27] 


1 


°4,3 


15KA^ 


60KAI 




°4,2 


35KAI 


UOKAl 


-175KAI 


°4,1 


+ K452A\ 


324KA\ + 1808-ff^ 


-405KAh - 2260KA\ 


°4,0 


U2KAI + 16KAI 


568KAI + 32KAI 


-710A A3 - 40K A\ 


°4,-l 


189KA 1 2 + W76KA\ 


756KAI + 4304KAI 


-945KA\ - 5380A~Ai 


Cl,-2 


215KAI 


860KAI 


-im$KA\ 


/~*l-t- 
°4,-3 


315KAI 


1260KAI 


-1575KA2 


°4,-4 


KA\ 


4KA\ 


-5KA\ 


°2,3 


-24KA\ 


96(2 - 9rj)Al 


24(1 - 39t;)^^ 


^1+ 

°2,2 


-4QKA\ 


8(4 - 123r7)A| 


-8(103 + 24377)^J 


°2,1 


-U4KAI - 256KA\ 


96(5 - Urj)Al - 64(13 + U4rj)Al 


48(3977 - U)Al + 64(42r7 - 167)Ai 


/~fl+ 
°2,0 


-Va2KA\ + 2KA\ 


16(34 - 123t])4J - 8(2 + Vrj)A\ 


8(69t7 - 65)Al + 2(39t? - 73) A 


r l+ 
°2,-l 


-504KAI + 224KA\ 


96(45?) - 22)^2 - 64(29 - 120?7)4i 


8(543 - 87377)^2 + 32(2377? - U5)A\ 


°2,-2 


-4KA\ + 32KAI 


8KA\ + 8(100 - 171?7)A3 


28KA\ + 336(128 + 2,ri)A\ 


r l+ 
°2,-3 





96(11 - 1677)^2 


96(23 - 1877)^2 


°0,3 


SOKAl 


24(29 - 15??)^ 


6(7577 - 121) A 


°0,2 


-SQKAl 


120(1 + 3r))Al 


-30(5 + $ri)A\ 


°0,1 


2MKA\ - 24KA\ 


24(2357? - 33)A^ - 96(37 - 35r?)Ai 


6(93 - 8237?)^ + 24(149 - 14377)Ai 


°0,0 


3KA{ - 6KA3 


-12KA{ + 24(13 - 4977)^3 + 12(3 - r))A\ 


9KA{ + 6(53 - 19977)A^ + 12(3 - 77)^5 



Table X: This table details the relevant coefficients for evaluating every nonzero Ck,j in Eq. (|41[) for h 1+ . It is worth mentioning 
that for h + every Sk,j are zero. 



Harmonic 


cos [37] 


cos [7] 


nlX 
J 4,3 


60KAI 


-607^ 


nlX 
J 4,2 


UOKAl 


-140^^ 


nlX 
J 4.1 


Z24KA\ + 1808A"Ai 


-324i^A + 1808^,41 


nlX 
J 4.() 


568KAI + 32KAI 


-568^^ + 327^ 


nlX 
J 2,3 


24KA\ 


24(7 - 6977)^2 


nlX 
J 2,2 


49KA\ 


-8(109 + 11777)^3- 


nix 
J 2,l 


~24KA\ + 5UKA\ 


24(3777 - 7)A\ - 32(585 + 69t7)Aj 


nlX 
J 2.() 


16KAI - 16KAI 


-112(1 + 977)^^ - 16(11 - 377)Aj 


nix 


-48KAI 


487^2 


nix 
°0,2 


-48KAI 


48KA\ 


nix 
J 0.1 


-487^^2 + 192KA\ 


48KA\ + \92KA\ 


nl X 
J -2,3 


-\2<dKA\ 


24(12177 - I3I7M2 


nl X 
°-2,2 


8KA\ - 2327^3- 


-40KA{ - 8(203 - 9377)^ 


nl x 
"-2,1 


672A"A 2 - 1696A:Ai 


24(49577 - 277)^ + 32(133 - 13577)Ai 


nl X 
°-4,4 


-4KA\ 


4KA\ 


nl X 
°-4,3 


-1260A'^2 


12fflKA\ 


nl X 
S -4,2 


-8&QKA\ 


860KAI 


nl X 
°-4,l 


-7b&KA\ - 4304A'Ai 


7hQKA\ + 4304A-Ai 



Table XI: This table details the relevant coefficients for evaluating every nonzero Skj in Eq. (|41[) for h lx . It is worth mentioning 
that for h lx every Ckj are zero. 



expressions for each coefficient become significantly more complex. It is because in addition to the harmonics of the 
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Harmonic 


cos [27] 


1 


°2,2 


-BAg 


-3BA& 


°2,1 


(1377 - 31)^ - (10 - 9rj)A& - 2(14 - 17r,)A? 2 


3(1377 - 31)A$ - 3(10 - 9ri)Au - 6(14 - 17t7)A^ 


n N+ 
°2,0 


2(577 - 3)^ + 8(577 - 2)^ 


6(577 - 3)Af + 24(5?7 - 2)A$ 


°2,-l 


(65t? - 51)Af - (32 - 4377)ylf 1 - 2(18 - 2377)^ 


3(65t7 - 5i)A? Q - 3(32 - 43r/)A^ - 6(18 - 2377)^ 


r ,N+ 
°2,-2 


-12KA? + (3377 - 32)A^ + (5377 - 48)A^ 


-36KA? + 3(3377 - 32) Ag + 3(53?? - 48)Af 


°2,-3 


(377 - 4)< 


3(377 - 4)Af 


°0,3 


2(2- 77)^ 


6(77 - 2)< 


s~i-N + 
°0,2 


(5 - ArftAg + 4(12-117?) Ag 


4(4?7 ~ 5)Af - 4(12 - llrj)A^ 


°0,1 


2(43 - 49r))A? + 2(41 - MrfA^ + 4(28 - 25t?) A? 2 


2(49?7 - 43)A^ - 2(41 - 3477)^ - 4(28 - 2577)Af 2 


°0,0 


12KAj + 2(21 - 1977)^ + 12(6 - 5-q)Ag 


-127\A^ - 2(21 - 19?7)Af - 12(6 - 5rj)Ag 


ciV+ 
°2,0 


-ABAg 


— I2BA5 


ciV+ 
°0,1 


2BA% 


-2BA% 


qN+ 

°-2.1 


2BA% 


6BA4 



Table XII: 1PN corrections to the coefficients for h N+ in Eq.(j39}. The table contains all coefficients to evaluate each nonzero 
Ck,j and Sk.j- Any other coefficients for these expressions are zero. 



Harmonic 


cos [7] 


r ,Nx 
°2,0 


16-BAf 


/"tiVx 
°2.-l 


8BA4 


oJVx 
J 2,2 


-ABA® 


qNx 
S 2,l 


4(1377 - 31)A? - 4(10 - 9t7)^S - 8(14 - 1777)Af 2 


qNx 
S 2,0 


8(577 - 3)Af - 32(2 - 577)^) 


qNX 

°-2,3 


4(4 - 377X 


qNX 

°-2.2 


4SKA? + 4(32 - 3377)^ + 4(48 - 5377)^ 


qNX 

°-2.1 


4(51 - 6577)Af + 4(32 - 437?)^ + 8(18 - 23rj)Ai 2 



Table XIII: 1PN corrections to the coefficients for h in Eg, (|39[) . The table contains all coefficients to evaluate each nonzero 
Sk,j and Ck,j- Any other coefficients for these expressions are zero. 



angle 7 there appear the harmonics of To, too. 

To keep our tables simple we introduce the following short-hand notations: 



c»+ = 


Cf c + 2 cos[2 7 ]H 








' ^feO 


q<Nx _ 








sr=^ cl cos[ 7 ] 




c u = 


C fe s 5 3 + sin[3 7 ]- 


, riO.5- 

r °feSl 


_ sin[7] , 


S° k 5+ = S° k g+sm[H 


+ 5° s 5 1 + sin[7] 


q,0.5x _ 


C°J 2 X sin [2 7 ] 






^• 5X =^ 2 Xsin[2 7 ] 




k + = 




r <l+ 


cos[2 7 ] + C+ , 


Sl + = Sl+ 4 cos [4 7 ] - 


r-Sj£ 2 cos [2 7 ] + S 


<%* = 


C^ 3 cos[3 7 ]H 


r^ilx 
~ °fcCl 


cos [7] , 


^ X =S% a cos [3 7 ] + 


Slci cos [7] ■ 



For open orbits we separate the higher harmonics in To instead of 7, and hence we write 



(46) 



C f = J2{o? f sm(ZTo) 



bf F cos (ZT„) 



S F = 



1=0 



\af F sin (IT ) + bf F cos (IT ) 



(47) 



where F is also a multi- index as for the elliptic case containing also the index of the harmonics of 7. Tables IXI VIIXXI Vl 



contain all the relevant a l F , b l F 



and b l F coefficients. Any other ai or bi coefficients are zero. 
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Harmonic 


cos 2To 


sin 2To 


1 


n N 
°3C2 


A-2 








^2C2 


A% 








^1C2 


5A? 







°0C2 


Af 





-A? 


r,N 

°30 


3< 








s~iN 
^20 


3^ 








/-iN 
°10 


15< 







Coo 







Af 


qN 
<->3C2 





-A? 





qN 
S 2C2 










qN 
°1C2 










qN 

d'M 





-ZA" 





qN 
"J 20 





-3Af 





C N 
"10 











Harmonic 


cos 2To 


sin 2To 


°3C1 





4A? 


<^2C1 





AA? 


rtN 





20A^ 


/-«JV 
°0C1 





4A? 


qJV 
D 3C1 


4A? 





oJV 
' 5 2C1 


4A? 





C JV 
"->1C1 


20A? 






Table XIV: This table contains the Newtonian value of all relevant a^ F , af F , bj F , and bf F coefficients needed for evaluating 
every nonzero Cp and Sp in Eq. (|43[) . The left table contains the coefficients for h , and in the right one the coefficients for 
h Nx can be found. 



For open orbits there are no difficulties when evaluating series expansion, in the final formulae there are two 
different types of terms: one contains higher harmonics of the true anomaly parameter, but there are also terms 
~ X cos kx/x sm kx- In Table IXXIVI the coefficients Ck x and Sk x of these terms are presented. 
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Harmonic 


COS [3To] 


cos [To] 


sin [3To] 


sin [To] 


n 0.5 + 
°6S3 


-4A?' 5 





A o. 5 





^0.5 + 
°5S3 







-A° 2 5 





^0.5 + 
°4S3 


— 12( J 4?" 6 + A"' 5 ) 


44 1 ' 5 


-11 A /" + 2A 3 h5 


11A ( 1 5 


^0.5 + 
°3S3 


-4(6AS>' 5 + < 5 ) 


12^' 5 


-5(2^' 5 - Al 5 ) 


23Ag s 


^0.5+ 
°2S3 


-12(A?- 6 +A§- 5 ) 


4(4y4?- 5 + 3,4§- 5 ) 


-13A?' 5 + I4A3' 5 


2(28^?' 5 + llAf 5 ) 


/-iO.S+ 
°1S3 


-12A° 2 5 


4(9A§' 5 + 


15A^ 5 


77 A% 5 + 13A° 4 5 


^0.5+ 
°0S3 




12« 5 + A% 5 ) 


7A o.s 


Ai' 5 + 42A% 5 


n 0.5 + 
°6S1 


-2(M?' 5 





5A ( l 5 





°5S1 


-60A§- 5 





-5A° 2 - 5 





^0.5 + 
°4S1 


-60(A?- 5 + ^' 5 ) 


2(L4? 5 


5(-lL4?- 5 +2A 3 '' 5 ) 


-25A?' 5 


^0.5+ 
°3S1 


-20(6^' 5 +^' 5 ) 




-25(2A§- 5 -ylj 5 ) 


-29A° 2 5 


/"i0.5+ 
°2S1 


-60(A?' 5 +^' 5 ) 


20(4A?- 5 + 3Ag' 5 ) 


-5(13A?- 5 - 14,4§- 5 ) 


2(-6^?' 5 + 7A§- 5 ) 


^0.5+ 
°1£T1 


-6(L4§- 5 


20(9 J 4g' 5 + ^- 5 ) 




81A° 2 5 + 17A° 4 - B 


°0S1 


-2(L4?- 5 


60(A?' 5 + A^ 5 ) 


7A o.s 


21Af B + 50A° 3 5 


C0.5+ 

°6S3 







4< 5 





C.0.5+ 

°5S3 







12A^ 5 





q0.5+ 

°4S3 


-11A?' 5 + 2A§' 5 


1L4?- 5 


12(A?' 5 +A% 5 ) 


-4A ( l 5 


q0.5+ 

°3S3 


-5(2A§' 5 - AT°) 


23yl§- 5 


A(6A'i- 5 + A° 4 5 ) 


-12A'i 5 


o0.5+ 
°2S3 


(-134?- 5 + 14A§- 5 ) 


2(5yl?- 5 + 1L4§' 5 ) 


12(A?- S + A§- 5 ) 


-4(2A?' 5 + 3A° 3 5 ) 


q0.5+ 
D 1S3 


VoAY° 


Aa 1 ' 5 + 13ASP 


12A§' 5 


-4(3y4 ( 2 ) ' 5 + A° 4 B ) 


q0.5+ 

°6S1 


5A?' S 





20 J 4°' 5 





q0.5+ 

°5S1 







6(L4|' 5 





q0.5+ 

°451 


-5(11A ( 1 ) ' 5 - 2A|P) 


-25A?' 5 


60(A ( 1 '' 5 + Ag' 5 ) 


-2(L4?' 5 


C0.5+ 

°3S1 


-25(2A 2 '' 5 -42' 5 ) 


-2QAY° 


20(6A§' 5 + Af 5 ) 




O0.5+ 

°2S1 


5(-KU?' 5 + 14A!p) 


2(-23A?' 5 + 7A§' 5 ) 


60(K 5 + A% 5 ) 


-20(2A° 1 5 + 3A° 3 5 ) 


Q 0.5+ 
°1S1 




-21A%- 5 + 17A° 4 5 


60AV 


-20(3A 2 )5 + Al 5 ) 



Table XV: This table contains the relevant af F , af F , bf F , and bf F coefficients for evaluating the explicit form of h°' 5+ in 
Eq. (|44)l . 
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Harmonic 


COS [3To] 


cos [To] 


sin [3To] 


sin [To] 


^0.5 X 
°6S2 


-4A?' 5 





-16.4?' 5 





^0.5X 
°5S2 







-48A^ 5 





^0.5 X 

°4S2 


4(11^?- 5 - 2A% 5 ) 


-4A , 5 


-48(A?- 6 + AD 


16A'l 5 


/-i0.5 X 
°3S2 


20(2yl^ 5 - A2' 5 ) 


-20A° 2 3 


-16(6A° 2 - 5 + A° 4 - 5 ) 


48A° 2 5 


^0.5 X 
°2S2 


4(13A?' 5 - I4A3' 5 ) 


-8(6 A f l 5 + 5A% 5 ) 


-48(A;' 5 + A% 5 ) 


16(4A't 5 + 3A$ 5 ) 


/i0.5x 
°1S2 


-60A^' 5 


-4(39A^ 5 + 7A° 4 5 ) 


-48A° 2 5 


16(9A% 5 + A° 4 5 ) 


^0.5 X 

°0S2 


-28A?' 5 


-4(3Af B + 22A°/°) 


-16.4? 5 


48(Af 5 + A*i- 5 ) 


rrO.SX 
J 6S2 


-16A?' 5 





4A , 5 





o0.5x 
J 5S2 


-48A^' 5 










crO.SX 
J 4S2 


-48(4?' 5 + A§' 5 ) 


IQAY 


4(-llA$- 5 +2AY j ) 


4A ( l 5 


o0.5x 
J 3S2 


-16(6A° 2 B + Af s ) 


48A° 2 5 


-2Q(2A° 2 5 - A° 4 S ) 


20A9,- 5 


o0.5x 
°2S2 


-48(A?- 5 + A° 3 - 5 ) 


16(2A?' 5 + 3A% 5 ) 


4(-13yl?' 5 + UA° 3 5 ) 


-8(A ( l 5 -5A$ B ) 


cO.Sx 
°1S2 


-48A° 2 5 


16(3A° 2 5 + Al 5 ) 


60A° 2 5 


4(-3A° 2 5 + 7A° 4 5 ) 



Table XVI: This table contains the relevant af F , af F , bf F , and bf F coefficients for evaluating the explicit form of h' 
Eq. (|44l) . 



Harmonic 


cos [4T ] 


cos [2To] 


1 


n i+ 

°7C4 


15KA\ 








°6C4 


■$hKA\ 








°5C4 


81KA\ + 452KA\ 


-2AKA\ 





r<i+ 

°4C4 


\A2KA\ + 8KA\ 


-40KAI 





°3C4 


189KA\ + W76KA\ 


~1UKA\ - 256KA\ 


30_ftr^2 


°2C4 


215KAI 


-152KAI+2KAI 


30KAI 


°1C4 


315KAI 


-504KAI + 22AKA\ 


2UKA\ - 24KA\ 


°0C4 


KA\ 


-4KA\ + 32KAI 


3KA\ - 6KAl 


°7C2 


mKAl 








°6C2 


UOKAl 








°5C2 


324KAI + 1808A:Ai 


96(2 - 977)^2 





r<l+ 
°4C2 


568KAI + 32KAI 


8(4 - 123r?)Aj 





r l+ 
°3C2 


756KAI + 4304A:Ai 


96(5 - 1477)^^ - 64(13 + lUr))Al 


24(29 - 1577)^2 


f-fl+ 
°2C2 


moKAl 


16(34 - 123rj)Al - 8(2 + 977)^5 


120(1 + 3^ 


r l+ 


1260KAI 


96(33 - &lf))A\ + 64(29 - 120*7)^4 


24(23577 - 33)A 2 - 96(37 - 3577)4i 


^.1+ 
°0C2 


4KA\ 


8KA\ +8(100- 171ri)A\ 


-12KA\ - 24(13 - 4977)^ - 12(3 - 77)^5 


°70 


-75KAI 








°60 


-VTbKAl 








r i+ 

°50 


4D5KAI + 2260KAI 


24(1 - 3977)^2 





WL+ 
°40 


-7WK Ak - 4QKA\ 


-8(103 + 2477)^4^ 





r l+ 
°30 


-945KA\ - 5380KAI 


48(3977 - 11) A\ - 64(167 - 4277)^! 


6(75?7 - 121)Al 


^1+ 
°20 


-1075KAI 


8(69?7 - 65),43 - 2(73 - 39r?)^ 


-30(5 + 9?7)Ai 


r l+ 


-1575KAI 


72(91 - 121??)A2 - 32(145 - 23777)4i 


6(93 - 82377)A^ + 24(149 - 143r7)A^ 


°00 


-5KA{ 


Z2KA\ + 8(128 + 3tj)^ 


9KA\ + 6(53 - 199t7)AJ + 12(3 - 77)^5 



Table XVII: This table contains all relevant coefficients for evaluating every Cf in Eq. (|45|l for h 1 



Harmonic 


sin [4T ] 


sin [2To] 


cl+ 

°7Ci 


-15KAI 





ql+ 


-3hKA\ 





ci+ 

D 5C4 


-SlKAj - A52KA\ 


24KAI 


ol + 
°4C4 


~\42KA\ -%KA\ 


AOKAl 


ol+ 
°3C4 


-ISQKAl - 1076^1 


U4KAI + 256KA\ 


ol+ 
°2C4 


-HhKA\ 


152KAI - 2KA\ 


ql+ 


-315KAI 


504K A\ - 224KA\ 


cl+ 

°7C2 


-60KAI 





cl+ 
°6C2 


-UOKAl 





ol+ 
°5C2 


-324KA\ - 1808^^1 


96(9t? - 2)A\ 


ci+ 

°4C2 


-56&K A\ - 32KA\ 


8(123?7 - 4)A\ 


ol+ 
°3C2 


-75QKAi, - 4304^^1 


96(14?? - 5)^2 + 64(13 + l\4r))A\ 


ol+ 
°2C2 


-860KAI 


16(1237? - 34)^3 + 8(2 + 9r))Al 


ol+ 
°1C2 


-1260KAI 


32(87?; - 33)^2 - 64(29 - 120r/)Al) 


$70 


ToKA\ 





ol+ 
J 60 


\7$KA\ 





ol+ 
D 50 


-4WaKA\ - 2260K A\ 


24(39?? - 1)A\ 


ol+ 
°40 


710KAI + 40KAI 


8(103 + 2477)^3- 


ql+ 


945KAI + 5380KAI 


48(11 - 39rj)Al + 64(167 - 43tq)A\ 


ol+ 
J 20 


1075A'Aj 


8(65 - 6977)^3 4- 2(73 - 397?)^ 


C-1+ 
°10 


1575KAI 


24(219?7 - 8V) A\ + (145 - 237?7)32Ai 



Table XVIII: This table contains all relevant coefficients for evaluating every Sf in Eq. (|45p for h 



Harmonic 


sin [4T ] 


sin [2T ] 


°7C3 


60A"^ 





v~*lX 

°6C3 


140A"Aj 





°5C3 


324KA\ + 1808^^1 


24KA\ 


°4C3 


568KAh + 32KAI 


4<dKA\ 


silX. 
°3C3 


756KAI + 4304KAI 


-24KA\ + 544KA\ 


v~flX 

°2C3 


860KAI 


-1QKA\ + 16KAI 


1C3 


1260KAI 


-576KA], + 16967^^1 


/llX 
°0C3 


4KA\ 


-8KA\ + 232KA\ 


^lx 
°7C1 


-WKA\ 





°6C1 


-14QKA\ 





z^lx 

°5C1 


-324KA\ - 1808^^1 


24(7 - 6977)42 


/llX 
°4C1 


-568KAI - 32KA\ 


-8(109 + 117>7)A3 


^lx 
^301 


-75QKAI - 4304KA\ 


24(3777 - 7)^2 - 32(385 + 6977)Ai 


WLX 

°2C1 


-860 K Al 


-112(1 + 9rj)Al - 16(11 - 3r))Al 


°1C1 


-1260^2 


192(51 - 7777)^2 - 32(133 - 135r7)A| 


v~flX 

°0C1 


-AKA\ 


40KA{ + 8(203 - 93r?)^ 



Table XIX: This table contains all relevant coefficients for evaluating every Cf in Eq. (|45H for h 1 



20 



Harmonic 


cos [4To] 


cos [2 To] 


1 


gl X 
°7C3 


&QKA\ 








nl X 
J 6C3 


UOKAl 








q1 X 
°5C3 


324KA\ + 1808KA\ 


247^2 





nl X 
°4C3 


568KAI + 32KA\ 


40KAI 





ol X 
°3C3 


756KAI + 4304KA\ 


-24KA\ + 544KA\ 


-4&KA\ 


nl X 
°2C3 


86QKAI 


-WKAl + 16KA\ 


-4&KA\ 


ol X 
°1C3 


12Q0KAI 


-81&KA\ + 1696KA\ 


-48KAI - \92KA\ 


ol X 
D 7G1 


-GOKAl 








ol X 

°6C1 


-UQKA\ 








ol X 
°5C1 


-324KA\ - 1808JsTAi 


24(7 - 69r t )Ai, 





ol X 
°4C1 


-568K A\ - 32K A\ 


-8(109 + 11777)^1 





q1 X 
°3C1 


-75QKAI - 4304JsTAi 


24(3777 - 7)^2 - 32(385 + 6977)Ai 


48KA\ 


ol X 
°2C1 


-860KAI 


-112(1 + 9t))Al - 16(11 - 3ri)Al 


48KA\ 


ol X 
°1C1 


12Q0KAI 


48(73 - 187tj)A^ - 32(133 - 135??)Aj 


487^^2 + 192KA\ 



Table XX: This table contains all relevant coefficients for evaluating every Sp in Eq. (|45[) for h lx . 



Harmonic 


cos [2To] 


sin [2T ] 


1 


r<N+ 

°4C2 










°3C2 


-(31 + 13r7)Af - 4(10 - 9ri)A& - 2(14 - 17 V )A? 2 





2(2-77)^ 


r>N+ 

°2C2 


2(-3 + 57 7 )ylf - WKAg 


-4BAs 


4(5 - 4r))Ag + 28(1 - n)A? 


°1C2 


{89rj - 83)Af - (48 - 557?)An - 2(22 - 31?7)Ai2 


-2BA$ 


2(43 - 49tj)Ai + 4(20 - H^A* + 4(28 - 25r))A? 2 


r-iv+ 

°0C2 


-12KA{ - (28 - l5f])Ag - (17 - 29??)^^ 





12KA\ + 3(9 - 5 V )Ag + (33 - 31ri)Ag 


n N+ 
°40 


-3BA? 








°30 


-3(31 + 1377)Af - 12(10 - 9?7)Afi - 6(14 - 177?)Af 2 





2(77 - 2)A% 


n N+ 
°20 


6(577 - 3)Ag - 30KAg 


-12BAg 


4(477 ~ *)*$ - 28(1 - 77X 


n N+ 
°10 


3(89t7 - 83)Af - 3(48 - 55rj)A^ - 6(22 - 31rj)A? 2 


—6BA4 


2(4977 - 43)Afo - 4(20 - Y7ri)A? x - 4(28 - 2hrf)A? 2 


n N+ 
°00 


-3&KA\ - 3(28 - l5t})Ag - 3(17 - 29^)^ 





-12K A\ - 3(9 - fnfiAg - (33 - 3l-q)A% 



Table XXI: 1PN coefficients for h N+ in Eq.(|43} from the 1PN description of the motion. This table contains only the 1PN 
corrections to the coefficients needed to evaluate the Cp quantities. 



Harmonic 


cos [2To] 


sin [2To] 


1 


qN+ 
°4C2 





BA» 





qN+ 
J 3C2 





(31 + 13rf)A N + 4(10 - 9?7)Af a + 2(14 - 17n)A? 2 





qN+ 
°2C2 


-4BA% 


2(3 - 5r))Ag + WKAg 





qN+ 
°1C2 


-2BA? 


(83 - 897?)Af + (48 - 55??)^ + 2(22 - 31rj)A? 2 


-2BA% 


qN+ 
°40 





3BAg 





°30 





3(31 + 13rj)A? + 12(10 - 977)^ + 6(14 - 17rj)A? 2 





qN+ 
°20 


-12BAg 


6(3 - 5rj)Ag + 30KAg 





qN+ 
°10 


-6BA? 


3(83 - 89r))A? + 3(48 - 55rj)A^ + 6(22 - 3l77)Af 2 


6BA% 



Table XXII: 1PN coefficients for h N+ in Eq.l|35]l from the 1PN description of the motion. This table contains only the 1PN 
corrections to the coefficients needed to evaluate the Sp quantities. 



27 



Harmonic 


cos [2To] 


sin [2T ] 







-4BAg 







-4(31 + 13v)Aw - 16(10 - 977)^^ - 8(14 - 1777),4f 2 




16BA$ 


8(577 - 3)Af - 40if J < 




8BA% 


4(8977 - 83)A? Q - 4(48 - 5577)Afi - 8(22 - 31ri)A? 2 







-16KA\ - 12(10 - 9r))Ag - 4(17 - 2377) A$ 


s? x 


-4BA5 





qNx 


-4(31 + 13r7)Af - 16(10 - 9^)^ii - 8(14 - 17r7)Af 2 





S? x 


8(577 - 3)4^ - 40KAg 


-16BAf 




4(89t7 - 83)Af - 4(48 - 55??)^ - 8(22 - 31r?)^ 


—SBAf 



Table XXIII: 1PN coefficients for h Nx in Eq.j43]) from the 1PN description of the motion. 



State 


+ 


X 


7-harmonic 


cos 


[27] 


1 


cos [7] 


To-harmonic 


cos [2To] 


sin [2T ] 


cos [2To] 


sin [2T ] 


cos [2To] 


sin [2T ] 


C-t x 





3A 6 " 





9A» 


-12^ 





C 2x 





48Ag 





144A^ 


-192Af 





Ci x 










45< 


-60Ag 





C 0x 





24Ag 





72Af 







Sa x 


3A 6 " 





9Ag 








12Ag 


S 2x 


4&A$ 





144A^ 








192A^ 


Si x 


ISA? 





45< 








60^ 



Table XXIV: This table contains the 1PN coefficients for both polarization states of Eg. (|43[) from the 1PN description of the 
motion which are multiplied by %■ The column under the sign " + " contains the coefficients for h N+ , and under the sign "x" 
the coefficients for h Nx can be found. 



